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CN ! Abstract 

Parallelograms are one of the basic building blocks in two-dimensional tiling. They have important 
applications in a wide variety of science and engineering fields, such as wireless communication 
networks, urban transportation, operations research, etc . Different from rectangles and squares, the coor- 
Oh . dinates of a random point in parallelograms are no longer independent. As a case study of parallelograms, 



the explicit probability density functions of the random Euclidean distances associated with rhombuses 
are given in this report, when both endpoints are randomly distributed in 1) the same rhombus, 2) two 
parallel rhombuses sharing a side, and 3) two rhombuses having a common diagonal, respectively. The 



J> , accuracy of the distance distribution functions is verified by simulation, and the correctness is validated 

■ by a recursion and a probabilistic sum. The first two statistical moments of the random distances, and 



the polynomial fit of the density functions are also given in this report for practical uses. 

Index Terms 

Random distances; distance distribution functions; rhombuses 

I. The Problem 

Define a "unit rhombus" as the rhombus with an acute angle of 9 — | and a side length of 
1. Picking two points uniformly at random from the interior of a unit rhombus, or between two 
adjacent unit rhombuses, the goal is to obtain the probabilistic density function (PDF) of the 
random distances between these two endpoints, as illustrated in Fig. \\\ 

There are four cases depending on the geometric locations of these two random endpoints, 
when rhombuses are adjacent and similarly oriented, as shown in Fig. [Q i.e., AB that are within 
the same rhombus; RS that are inside two parallel rhombuses sharing a side; PQ and MN that 



Fig. 1: Random Points and Distances Associated with Rhombuses. 



are inside two rhombuses sharing a common diagonal. Here PQ and MN are two different 
cases, and in the following, we refer to them as long (long-diag) and short diagonal (short-diag), 
respectively. The next section gives the explicit PDFs for the above four cases. 

II. Distance Distributions Associated with Rhombuses 

A. \AB\: Distance Distribution within a Rhombus 

The probability density function of the random Euclidean distances between two uniformly 
distributed points that are both inside the same unit rhombus is 



f Dl (d) =2d{ 



* (l + f ) -n- 1 # - 




9^ 



2tt 
9^ 



< d < & 
f<d<l 
l<d<V3 








The corresponding cumulative distribution function (CDF) is 
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B. \RS\: Distance Distribution between Two Parallel Rhombuses Sharing a Side 

The probability density function of the random distances between two uniformly distributed 
points, one in each of the two adjacent unit rhombuses that are parallel to each other, is 
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The corresponding CDF is 

-± (4* + d±) sin- f d + ( - I) + |rf3 + ^ 
- 94d2+15 -V4d^3 



72 

2d 4 • -1 
3^3 Sm 2d 



2d 2 _ d 4 \ • -1 , f 2d 2 , tf 4 \ • -1 V3 i ( 1 j4 



16. 

9 



55 

48 



"q _ 23 

48 



< d < & 



f<d<l 



l<d<V3 



V3<d<2 



2<d<^7 

otherwise 
(4) 



C. \PQ\: Distance Distribution between Two Long-Diag Adjacent Rhombuses 

The probability density function of the random distances between two uniformly distributed 
points, one in each of the two adjacent unit rhombuses that have a common long diagonal, is 
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The corresponding CDF is 
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ZX \MN\: Distance Distribution between Two Short-Diag Adjacent Rhombuses 

The probability density function of the random distances between two uniformly distributed 
points, one in each of the two adjacent unit rhombuses that have a common short diagonal, is 
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Note that although unit rhombuses are assumed throughout (d)-©, the distance distribution 
functions can be easily scaled by a nonzero scalar, for rhombuses of arbitrary side length. For 
example, let the side length of such rhombuses be s > 0, then 

F sD (d) = P(sD <d)= P(D <-) = F D {-). 
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Therefore, 
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III. Verification and Validation 



A. Verification by Simulation 

Figure |2] plots the probability density functions, as given in (QQ), ©, © and ©, respectively, 
of the four random distance cases shown in Fig. CD Figure |3] shows a comparison between the 
cumulative distribution functions (CDFs) of the random distances, and the simulation results by 
generating 1, 000 pairs of random points with the corresponding geometric locations as illustrated 
in Fig. [U Figure |3] demonstrates that our distance distribution functions are very accurate when 
compared with the simulation results. 
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Fig. 2: Distributions of Random Distances Associated with Rhombuses. 
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Fig. 3: Distribution and Simulation Results For Random Distances Associated with Rhombuses. 



B. Validation by Recursion 

When looking at Fig. \T\ we find that the four adjacent rhombuses together resemble a larger 
rhombus, with a side length of 2. According to ©, the distance distribution in the large rhombus 
is f2D(d) = hfDiif)' On the other hand, if we look at the two random endpoints of a given 
link inside the large rhombus, they will fall into one of the four individual cases: both endpoints 
are inside the same small rhombus, with probability |; the two endpoints fall into two parallel 



s 



rhombuses, with probability |; and the two endpoints fall into two diagonal rhombuses (either 
long or short-diag), both with probability |. Thus the distance density function for the large 
rhombus can be given by a probabilistic sum, f 2 D{d) = ifD 1 (d) + ^f Dp (d) + lfD LD (d) + lf DsD (d), 
where foM), fD P (d), fD w (d) an d fD SD (d) are given in (0Q), ©, © and ©, respectively. 
To confirm that the above two definitions of f2D{d) are equivalent, i.e., \fo l {d) + \fo P {d) + 
lfD hU (d) + §/i) SD (d) is equal to ^/di(§)> we verify them mathematically as follows. 
1) < d < ^: 



\fm(d) = \ 



- H f= d 2 d + — 

3 9VSJ 3 ^ 



:/dp(^) = d 



3 \3 9V3J 



Thus, 



d 
2 



4 2it \ fd 

+ 



3 9^7 V 2 



16 /d 

T 1 2 



2tt 

71 



2) ^ < d < 1: 



8 / d 2 \ . ,V3 /4 10vr\ 2 16, 10 / — — - 2vr 

— 1 H sin -1 — + = d 2 d H V4d 2 - 3 = 

v/3 V 3 y 2d \3 9V3J 3 3 ^3 



^fD P {d) = d 



•-^=(d 2 + 2) sin 



, V3 f 8vr 2\ 4 , 11 5 2tt 

1 ^— + — = - - d 2 + -d V4d 2 - 3 + -= 

2d V9V/3 3^ 3 6 y/3 



1 



^ / 1 7^ ^ J 2 1 r /A 

= d > q/d S dW 



4 V3 9^/ 



d 
4 



8d 2 . , >/3 / 1 IOtt \ l2 2 

— - sin^ 1 — — + I d 2 + -V4d 2 - 3 

3v/3 2d V3 9V3J 3 



9 



Thus, 



hD{d) = \f Dl (d) + lfD P (d) + lfD LD (d) + lfD s M 



d 

2 



4 2tt \ fd\ 16 fd\ 2tt 



3 9V3J V 2 / 3 V 2 / \/3 



\ fD A ] - 



3) l<d<VS: 



\fnM) = \ 



t V3 (2 2tt \ 



4 / d 2 \ . V3 /2 

— ^= 1 sin — - — - 

V3 V 3 ; 2c? V3 



2tt 



_ d 2 + V4d 2 - 3 - ^= - 1 
9y/3J 3^ 



2/Dp(rf) = d 



4d 2 ,V3 2 2tt \ l2 8 , V4rf 2 - 3 2vr 1 

— - sin -1 — — ■ + I — ) d 2 - -d H 1 = + - 

3v/3 2d \3 9^3/ 3 3 3^ 2 



4d 2 . V3 ( 7T 

3^3 Sm 2d \3>/3 



1 d 



2i ? 



/ 4d 2 8 \ . ,V3 (I 2- 



8tt 



2, ' V3'9V3/ 2 + 3 rf - 3Vi ^ + i7I +1 



Thus, 



d 
2 



I 2- \ /dV 16 /d\ 2 k 



3 9^3/ V 2 / 3 V 2 / y/3 
4) V3<d<2: f Dl (d) = 0, and 



2 -M 2 , 



j/op(d) = d 



2 d 2 \ , VS ( 2 d 2 \ ,V3 A 



7 



+v / d 2 ^3 + — V4d 2 -3 + ^ 

12 4 3^ 



2d ' Vv^ ' 3/37 ' "' V3 Ov^y ^ 3^ 

3 2tt 



7T 



^/d ld (^) = \ 



7T 



4 /d 2 \ ,V3 fl 

— 2 sin" 1 — — — h I — 

VS\3 J 2d \3 9VS 



10 



g/D SD (cO - ^ 



8 . i y/3 ,o 8 , 

—= sin — - a H — a 

y/3 d 3 



-Vd 2 - 3 = - 4 

3 3VS 



Thus, 



d 

2 



73 V 3 J 8111 2(d/2) + \3 9VS 



10 / (d 

+ T\ 4 2 



3-^ 
\/3 



2 fD A ] - 



16 /d 
T \2 



5) 2 < d < V7: f Dl (d) = f DsD (d) = 0, and 



^fD P {d) 



d 



2 d 2 \ ( . _! v/3 , . _! >/3 

— p= p= sm — - + sm — - 

y/3 3V3J 1 2d d 



4 + ^ 1 ^V(^- i l rf 2 + -v / 4d^3 
2d d J \Qy/3 3 J 12 



v/d^3 2tt 



3^ 4 



gfD LD (d) 



d 
4 



7T 



4 ^ 2 \ . iV3 2d 2 . _,VS . 
— 2 sin^ 1 -—- H = sm 1 — + 1 

y/3 V 3 / 2d 3^ d V 3^ 



2 m 4vr 



+ 3 



7 



d 2 - -V4d 2 - 3 



Thus, 



d 

2 



4 / (d/2f 
V 3 



sm 



d 



2tt 
9^ 



\ fD A ] - 



+ \ 4 



-3-^-1 



3^ 



11 
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In summary, we have f2o(d) = |/di(|) by recursion, and the probabilistic sum hf^d) + 
\}D P {d) + |/d ld (^) + |/d S d(^) * s ec l ua l t0 |/di(|) i n ai l tn e cases discussed above. The results 
are a strong validation of the correctness of the distance distributions that we have derived. 

IV. Practical Results 

A. Statistical Moments of Random Distances 

The distance distribution functions given in Section [XT] can conveniently lead to all the statistical 
moments of the random distances associated with rhombuses. Given fD x {d) in (0Q), for example, 
the first moment (mean) of d, i.e., the average distance within a single rhombus, is 
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and the second raw moment is 
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from which the variance (the second central moment) can be derived as 
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0.0708017742. 



When the side length of a rhombus is scaled by s, the corresponding first two statistical 
moments given above then become 
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TABLE I: Moments and Variance — Numerical vs Simulation Results 



Endpoint Geometry 


rDrvbim 


M D 


M D 


Varjy 


Within a 
Single Rhombus 




0.5123783359s 


0.3333333333s 


0.0708017742s 2 


Sim 


0.5137344650s 


0.3356749448s 


0.0717518443s 2 


Between two 
Parallel Rhombuses 


fD P (d) 


1.0750863337s 


1.3331823503s 


0.1773717254s 2 


Sim 


1.0749140141s 


1.3318514164s 


0.1764112787s 2 


Between two Long-Diag 
Adjacent Rhombuses 


fD hu (d) 


1.7570796617s 


3.3330145688s 


0.2456856311s 2 


Sim 


1.7575618714s 


3.3319543736s 


0.2429306419s 2 


Between two Short-Diag 
Adjacent Rhombuses 


/d S d (d) 


1.1150961004s 


1.3332064091s 


0.0897670959s 2 


Sim 


1.1156689048s 


1.3341275970s 


0.0894104918s 2 



TABLE II: Coefficients of the Polynomial Fit and the Norm of Residuals (NR) 



PDF 


Polynomial Coefficients 


NR 




10 s x [-0.000166 0.002955 - 0.024296 0.122501 - 0.423700 
1.065297 - 2.013141 2.916117 - 3.273038 2.858823 - 1.941184 
1.018567 - 0.408430 0.123035 - 0.027166 0.004245 - 0.000446 
0.000029 - 0.00000116 0.0000000892 - 0.0000000000648] 


0.095901 


fD P (d) 


10 5 x [0.000019 - 0.000513 0.006194 - 0.045754 0.231329 
-0.847940 2.328313 -4.879792 7.881444 - 9.836044 9.45481 
-6.93816 3.828928 - 1.553959 0.4491195 - 0.088213 0.010963 

-0.000782 0.0000499 - 0.0000000863 - 0.0000000025988] 


0.059485 




10 4 x [0.0000002 - 0.0000074 0.0001242 - 0.0012723 0.0089355 
-0.045576 0.174538 - 0.511656 1.160176 - 2.042559 2.78728 
-2.927992 2.339236 - 1.395868 0.606494 - 0.185173 0.037762 
-0.004740 0.0003293 - 0.0000098938 0.0000000733372] 


0.011340 


fD 3D (d) 


10 Y x [-0.000022 0.000493 - 0.005116 0.032591 - 0.142436 
0.452238 - 1.077816 1.965337 - 2.77016 3.029748 - 2.567377 
1.674771 - 0.831348 0.308601 - 0.083558 0.015932 - 0.002034 
0.000161 - 0.00000698 0.00000013178 - 0.0000000006098] 


0.147836 



Table U lists the first two moments, and the variance of the random distances in the four cases 
given in Section UH and the corresponding simulation results for verification purposes. 

B. Polynomial Fits of Random Distances 

Table HU lists the coefficients of the degree-20 polynomial fits of the original PDFs given in 
Section UIl from d 20 to d°, and the corresponding norm of residuals. Figure @](a)-(d) plot the 
polynomials listed in Table |n] with the original PDFs. From the figure, it can be seen that all 
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(a) Within a Single Rhombus 



(b) Between two Parallel Adjacent Rhombuses 





(c) Between two Long-Diag Adjacent Rhombuses 



(d) Between two Short-Diag Adjacent Rhombuses 



Fig. 4: Polynomial Fit of the Distance Distribution Functions Associated with Rhombuses. 



the polynomials match closely with the original PDFs. These high-order polynomials facilitate 
further manipulations of the distance distribution functions, with a high accuracy. 

V. Conclusions 

In this report, we gave the closed-form probability density functions of the random distances 
associated with rhombuses. The correctness of the obtained results has been validated by a 
recursion and a probabilistic sum, in addition to simulation. The first two statistical moments, 
and the polynomial fits of the density functions are also given for practical uses. 
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